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Roughly observed, some random phenomena seem scale invariant. This is the case 
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1. Introduction 
t> 

O 
00 
q ! for the velocity field of turbulent flows or the (logarithm of) evolution in time of the 

price of a financial asset. However, a more precise empirical study of these phenom- 
ena displays in fact a weakened form of scale invariance commonly called multifractal 
scale invariance or intermittency (the exponent which governs the power law scaling 
of the process or field is no longer linear). An important question is therefore to 
construct intermittent random fields which exhibit the observed characteristics. 

Following the work of Kolmogorov and Obukhov ([9], [12]) on the energy dissi- 
pation in turbulent flows, Mandelbrot introduced in [TD] a " limit-lognormal" model 
to describe turbulent dissipation or the volatility of a financial asset. This model 
was rigorously defined and studied in a mathematical framework by Kahane in [8]; 
more precisely, Kahane constructed a random measure called Gaussian multiplica- 
tive chaos. A natural extension of this work is to use Gaussian multiplicative chaos 
to construct a field (or a process in the financial case) which describes the whole 
phenomenon: the velocity field in turbulent flows (the price of an asset on a financial 
market). This extension was first performed by Mandelbrot himself who proposed 
to model the price of a financial asset with a time changed Brownian motion, the 
time change being random and independant of the Brownian motion. In [2J, the 
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2 HYDRODYNAMIC TURBULENCE AND INTERMITTENT RANDOM FIELDS 

authors proposed for the time change to take the primitive of multiplicative chaos: 
this gives the so called multifractal random walk model (MRW) (Bacry and Muzy 
later generalized the construction of the MRW model in |3j). The obtained process 
accounts for many observed properties of financial assets. 

The drawback of the above construction and of the MRW model is that the laws of 
the increments are symmetrical. In the case of finance, this is in contradiction with 
the skewness property observed for certain asset prices. In the case of turbulence, 
the laws of the increments must be nonsymmetrical: it is a theoretical necessity and 
stems from the dissipation of the kinetic energy ([7]). In light of these observations, 
we are naturally led to construct random fields which generalize to any dimension 
such process and which present multifractal scale invariance as well as nonsymmet- 
rical increments. 

We will answer a very natural question: how can one obtain a family of multi- 
fractal fields with nonsymmetrical increments by perturbing a given scale invariant 
Gaussian random field on R d ? Finally, in the last part we will mention the difficulties 
which arise in trying to construct an incompressible multifractal velocity field that 
verifies the 4/5-law of Kolmogorov with positive dissipation. 



2. Notations and preliminary results 

2.1. The underlying Gaussian field. Let dWo(x) denote the Gaussian white 
noise on M. d and ip : M. d — > [0, 1] denote a C°°, radially symmetric function worth 1 
for \x\ ^ 1 and for \x\ > 2. We also introduce a fixed correlation scale R > and 
a a number which satisfies 

d/2 < a < d/2 + 1. (2.1) 

We define the Gaussian field X g by the following formula: 

X g (x) = [ <p R (x - y) X ~l a+1 dW (y), (2.2) 

jR d f — y\ 

where we set the following notation: 

Using Kolmogorov's continuity criterion (see the standard book [13]), it is easy 
to show that (12. 2p defines a homogeneous, isotropic gaussian field which is almost 
surely Holderian of order < a — d/2. Note that condition (12. II) implies that the 
integrand in (12.21) is square integrable and the R d / 2 - a factor ensures that the field 
is dimensionless. 

Scaling property. Let e be a unitary vector and A > 0. We have the following identity 
in law: 



( Pn(y)y _ <PR(y-^e) (y - Ae). 
(law) J R d y \y\ d - a+1 \y-Xe\ 



X g (x + Xe)-X g (x) = s / (-^— - x )dW Q (y). 
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From the Gaussianity of the above law, we deduce that for all q > 0, there exists 

E(\X g (x + \e) - X g (x)\«) = alc g , 



c q > such that: 



with 

2 _ f ( l fR(y)y <Pr{v- Ae)(y- Ae) 2 

° Xe ~L [ \y\ d - a+1 ~ ly-Ael*-* )«v 

_ X 2a-d f , <Pr(*z)z _ <p R (\z - \e)(z - e) 2 
A" M d_Qf+1 |;z - e| d - Q+1 J 

^ (_A 2Q! - d / I z z — e 2 

xZo y R' J Rd ' \ z \ d - a+1 " \z - e\ d - a+1 ' 

We thus derive the following scaling 

A 



E(\X g (x + Ae) - X g (x)\*) ~ (^) 9(a - d/2) C g , 



where the constant C g is independent of e. One says that (X g (x)) x£R d is at small 
scales monofractal with scaling exponent a — d/2. 

A homogeneous and isotropic field (X(x)) xeR d is multifractal if there exists a non 
linear function ( q such that: 

E(\X{x + \e)-X(x)\«) ~ (h Cq C q . 

A— >0 it 

We call C q the structural function of the field (X(x)) xGVi d. 

2.2. Outline of the construction of multifractal vector fields from the field 

X g . Our construction is inspired by the work of Kahane in [S]. Let e > and X e (y) 
be a regular family of scalar Gaussian fields (not necessarily independent of dW ) . 
We consider a family of fields X e (with scalar components Xj! in the canonical basis) 
defined by: 

KM = f M* - v h x ~J a+ y Hy) - c ' dw o(y) ( 2 - 3 ) 

V \x-y\i + 

(\x — y\ e is defined in the next subsection and is given by a standard convolution). 
For an appropriate family X e , we show that it is possible to find constants C e such 
that X € tends to a non trivial field X (with scalar components X^ in the canonical 
basis) as e tends to 0. If one chooses X e independent of dW , we will see that this 
leads to a field X that extends the model introduced by Bacry in [2j and that 
has symmetrical increments. Thus, to obtain nonsymmetrical increments, we must 
introduce correlation between X e and dW . 

2.3. Notations and construction of the family X e . Let k R be the function 

■ T-tn? for Ixl ^ R, 
k R ( x ) = { W d/2 II-' 

otherwise. 



4 HYDRODYNAMIC TURBULENCE AND INTERMITTENT RANDOM FIELDS 

Let 9{x) be a C°°, non negative and radialy symmetrical function with compact 
support in |x| ^ 1 such that 



9(x)dx = 1. 

We define 9 e = -^O(-) and the corresponding convolutions: 

k* = 8 e *k R , |.| e = e *|.|. 

Let 7 be a strictly positive parameter and dW be a gaussian white noise on M. d . We 
consider the following gaussian field: 

X\y) = i [ kf{y-a)dW{a). 
Its correlation kernel is given by: 

E{X\x)X\y)) = ^p e/R (^), 

where p = k 1 * k 1 and p e = 8 e * 8 e * p. One can prove the following expansion 

p(x) =u d \n + — - + 4>(x), 
\x\ 

where uj^ denotes the surface of the unit sphere in M. d and is a continuous function 
that vanishes for \x\ ^ 2. We will note |.|* = inf(l, |.|) and, with this definition, the 
previous expansion is equivalent to: 

p <t>(x) 
e P(x) = 






One can also prove the following expansions with respect to e for e < R: 

*?(°) = % ( 2 - 4 ) 

with C = Jj u , < 1 rWi du and there exists a constant C\ such that 

p e/R (0) = to d lnj + C 1 + o(e). (2.5) 

In the sequel, we will consider the case 

<ydW = 7o(e)dW / + lidW u 

where dW\ is a white noise independant of dW and 70(e) is a function of e that will 
be defined later. Note that the integral in formula (12.31) has a meaning since dWo 
can be viewed as a random distribution. 



HYDRODYNAMIC TURBULENCE AND INTERMITTENT RANDOM FIELDS 



2.4. Preliminary technical results. We remind the following integration by parts 
formula for gaussian vectors (cf. lemma 1.2.1 in |llj): 

Lemma 2.1. Let (g, g±, . . . , g n ) be a centered gaussian vector and G : W 1 — > K a C 1 
function such that its partial derivatives have at most exponential growth. Then we 
have: 

n dG 

E(gG( 9l , ..., g n )) = ^ E(g gi )E( — ( 9l , . . . , g n )). (2.6) 

4 = 1 * 

From the above formula, one can easily deduce by induction the following lemma 
which will be frequently used in the sequel: 

Lemma 2.2. Let I G iV* be some positive integer and (g,gi, ■ ■ ■ ,g2i) Q> centered 
gaussian vector. Then: 

i 



E( gi ...g 2l e 9 ) = (J2 S ^ e 



fc=0 

where 



Sk,i= Yl ^2 E (99h)---E(gg i2k )E(g i2k+l g i2k+2 )...E(g i2l _ 1 g i2l ), 

{il,...,i 2 k}C{l,...,2l} 

where the second sum is taken over all partitions of {1, ... , 2l}\{ii, . . . , i 2k } in 
subsets of two elements {i2 P +i, ^+2}- 
Similarly, we get the following formula: 



i 



E(gi . . . g2i+ie 9 ) = (£ S ki i)e 



fc=0 

where 



Sk,i Yl S E (99h) ■ ■ ■ E (99i2 k +i) E (9i2k + 29i2k +3 ) ■ ■ ■ E (9 l2 i9i 2 i+i), 

{i 1 ,...,i 2k+ i}C{l,...,2l+l} 

Remark 2.3. In S k ,i (S k ,i), the summation is made of 2kl2 i-k {l _ k)l ( \ 2 k+iy^-Hi-kyJ 
terms, number we will denote by a k) i (a k j). 

We will also use the following lemma essentially due to Kahane ([8]). 

Lemma 2.4. Let (T, d) be a metric space and a a finite positive measure on T 
equiped with the borelian a -field induced by d. 

Let q : T x T :—>■ M. + a symmetric application and m a positive integer. Then we 
have the following inequalities: 



f e ^<i<*<a»«^'**W(ti)...d(T(t2 m ) < <x(T)(sup / e mq{t ' s) da{t)) 

JT 2m sGT Jt 



m 9(M).W-m2m-l 



(2.7) 
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e 



J'2m+l 



Ei ^ J<fe < 2m+ i rtj^daih) . . . da(t 2m +i) 



,s JT 



^ a{T) sup( / e q ^da{t)){ / e mq{s > t] e q{s > l) dait)) 2 ™- 1 . (2.1 



Proof. The proof of (\2.7\i can be found in [8]. Thus we just prove how to derive 
inequality ( 12. 81) from ( 12. 7ft . By integrating with respect to the first 2m variables and 
applying (12. 7p with the measure e g ^' t2m+1 ^da(t), we get: 



/ e Ei ^ i<fe ^ 2m+ i iVi^do-itx) . . . do-(t 2m +i) 

Jj^2m + 1 

2 n i 



= f do-(t 2m+1 ) f e^^<^^ q ^' tk) T[e q ^^" +l) da(t l )...da(t 2m ) 

JT JT 2m =1 

^ / da{t 2m +i){ [ e 9( *'* 2m+l) da (£)) (sup / e m « (8 '* ) e ff(t '* ara+l) dcr(t)) 2m ^ 1 

1277] Jr Jt a Jt 

^ a(T) sup( / e«6'W(t))( / e^'V^W^)) 2 ™" 1 . 

3. Construction of a four parameter family of multifractal, 
homogeneous, isotropic vector fields with non symmetrical 

increments 

In this section, we will suppose that d/2 < a < (d/2 + 1) A d and uJdll < d. 

We consider the field X e defined by formula ( 12.31) with 

X*(y) = l0 (e)X^y) +ll X{(y), 



where 

We set also 
and 



X-(v)= [ k?(y-<r)dWi(.<r), i = 0,l. 

JR d 
a = ((70(6)) 2 + 7l 2 )PeM0). 

7o (e) = 7o*( 



e , d -">dij 



R' 

Therefore, we introduce a slight correlation between X e and dWo (70(e) tends to 
as e goes to 0). 

3.1. Multiplicative chaos in dimension d. Multiplicative chaos or the "limit- 
lognormal" model introduced by Mandelbrot is a generalization of the exponential 
of a gaussian process. As mentioned in the introduction, it was defined rigorously by 
Kahane in [8] . The construction of Kahane was based on the theory of martingales 
and thus the generalized correlation kernel (here p(t — s)) had to verify a condition 
hard to verify practically (the a-positivity condition). Our construction is based on 
L 2 -theory and can be carried out without this condition. 
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We will construct the multiplicative chaos associated to the generalized correlation 
kernel p{^-jf-) defined in 2.3 and to some (positive) intermittency parameter 71 such 
that •yfujd < d. 

Let e be a positive number. Let B(M. d ) denote the standard borelian a-field; we 
want to consider the limit as e goes to of the random measures Q e ' 71 defined by: 

Q^(dy) = e -n*f(i/)-#s((*f&)) 3 ) dy 

This leads us to state the following proposition: 

Proposition 3.1 (Multiplicative chaos of order 71). There exists a positive random 
measure Q 11 (dy) independent of the regularizing function 9 such that: 

(1) for all A bounded in B(R d ), E(Q^(A)) = \A\. 

(2) Q 11 has almost surely no atoms. 

(3) Almost surely, Q 11 is singular with respect to the Lebesgue measure on all set 
A (with positive measure). 

If q is some positive integer and f : M. d — » R a deterministic function that satisfies 
the following condition: 



[ |/(yi)|-..|/WI n „ !^ dyx...dy 2q < 00, (3.2) 

then we have the following convergence: 



f(y)Q e ™(dy)^ Q J f(u)Q- A (<hn. 



We also have the following expression for the moments of J Rd f (y)Q 11 (dy) : 

VA;^2g, E ( f(y)Q^(dy)) k ) = f( Vl ) . . . f(y k ) f] d Vl ...dy k . 

(3-3) 
Moreover, the above formula ^3.3\) extends straightforwardly to the case of two func- 
tions f,g and two intermittency parameters 71,72 giving: 



x 



E[{ f(y)Q^(dy)) k ( g(y)Q»{dy)) 1 

= / f{yi)---f{yk)g{yk+i).-.g{yk+i) f\ .,.,, 

J ( R > + l^i<j^k |— R-|* 

11 \vi-v mw 11 ■tt-wiTfr., '**' 1 ••• dy *+'- 

l^i^k,j>k< Ft, I* fc+1 < i<j s? k+l \ J^ - I* 



We will call Q" 11 (dy) multiplicative chaos of order ^1 . 



HYDRODYNAMIC TURBULENCE AND INTERMITTENT RANDOM FIELDS 



Proof. We first start by considering a positive integer q and a function / that satisfies 
the corresponding integrability condition (13. 2p . Let e, e 1 be two positive numbers. By 
using Fubini, we get for all j ^ 2q: 



E(( f(y)Q em (dy)) j ( f(y)Q e '^(dy)) 2 ^) 



-i7?P,/«(0)-^T?P e7 H(0) / f( yi )...f(y 2q) X 

(Rd)2q 
e H*(Etl^(W)+a +1 ^(i«)) a ) dyi . . . d y 2q 

2A,f y i- y 3 



f e 7f<M— R^) 

f(yi) ■ ■ ■ fiihq) [ [ — — 1- 4 • • • ^ 



since p € /r(x) — > p(x) as e goes to 0. 
From this, we deduce that: 



E(([ f(y)Q £ni (dy)~ I f(y)Q e>ni (dy)) 

JR d JR d 






and therefore that J Rd f(y)Q eai (dy) is a Cauchy sequence in L 2<? that converges to 
some random variable Q 11 (f). For k ^ 2q, the moment E((Q J1 (f)) k ) is the limit as 
e goes to of E((Q e ^ 1 (f)) k ); from this one can deduce that the moments of Q' yi (f) 
are given by formula (13. 3p . 

For any bounded set A in B{R d ), consider / = 1^ and q = 1. Since 7^0^ < d, 
we deduce from lemma [231 that the integrability condition (13.21) is satisfied. Thus it 
follows from the proof above that Q eni (A) converges in L? to some random variable 
Q 11 (A). This defines a family of random variables (indexed by the bounded Borelian 
sets) that satisfies the following properties: 

(1) For all disjoint and bounded sets Ai, A 2 in B(R d ), 

Q 71 (A x U A 2 ) = Q 71 (A^ + Q 71 (A 2 ) a.s. 

(2) For any bounded sequence (A n ) n j> 1 decreasing to 0: 

Q J1 (A n ) -* a.s. 



n^oo 



By theorem 6.1. VI. in [5], there exists a random measure Q 11 such that for all 
bounded A in B(M. d ) we have: 

Q^(A) = Q^(A) a.s. 

Finally, one can easily show that the limit random variable Q 71 (/) is almost surely 
equal to J Rd f(y)Q" fl (dy). □ 

3.2. Convergence of X e towards a field X. 

In the sequel, (ej)j will denote the canonical basis (whereas (e J )j denotes 
the components of a vector e). In this subsection, we will prove the following 
proposition: 
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Proposition 3.2. Let a be such that d/2 < a < {d/2 + 1) A d and '-fi such that 
2 r )\ujd < a — d/2. There exists a field (X(x)) x€ ^d such that for all k and x±, . . . ,x k € 
M. d the following convergence in law holds: 

(X e ( Xl ), ..., X e (x k )) =► (X( Xl ), ..., X(x k )). (3.4) 

Let I be an integer such that one of the following conditions hold: 

(1) I is even and l^\u] d < a — d/2. 

(2) I is odd and (I + l)j((jj d < a — d/2. 

Let F R be defined by F R (y) = yj?(y) i >d- a +i ■ Then there exists C such that, for all x 
in IR d ; the random variables X^(x) have a moment of order 21 given by the following 
expression: 

E((Xi(x)f) = J2 a k ,C 2k f F R ( yi ) . . . F R {y 2k ){F R {y 2k+1 )f . . . {F R {y k+l )f 

fc=0 J(K d ) k + l 

11 1 vi-yj \il"d 11 wi-y, \ii\ud 11 \ yi-yj_ \^l^d yi '" Vk+l ' 

1 < i<j ^ 2fc I fi I* 1 < i < 2k I fi \* 2fc+l < i<j < fc+i | /j I* 

j>2fe 

(3.5) 
£((*'■(* + h) - Xi(x)f) =j^a Kl C 2k [ (F R ( yi ) - F R ( Vl -h))... 



k=0 

(Fi(y 2k ) - F R {y 2k - h)){F R {y 2k+1 ) - F R (y 2k+1 - h)f . . . (F R (y k+l ) - F R {y k+l - h)f 
11 \y.-y,al^ 11 \v,-vuH^ 11 ,*-„, • ,***> ' ' ' ^ (3 ' 6) 

1 < i<j ^ 2fc I /j I* 1 < i < 2fc I # I* 2fc+l < i<j < k+l \ ^ I* 

i>2fe 

Proof. Let 71 be such that 2^\ud < a — d/2. We set 

and define two auxiliary fields y e , Z e by the following expressions: 

y e (x) = f M* - 1/) , x ~l a+1 Q e ' yi (dy) (3-7) 

and 

Z e (x) = / <p R (x - y) . X ~l a+1 e-* x M- c *dW Q (y). (3.8) 

jR d f ~~ y\ 

Note that Z e (x) exists since X\ and dWo are independent with: 

JM d F y 
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We can compute, for all x in M. d , E(\X e (x) — Cy e (x) — Z e (x)\ 2 ) (cf. the more 
complicated computations in the proof of proposition 13.71) and derive the following 
limit: 

x t {x) - (cy e (x) + z e ( x )) * o. 

Thus, we must show that the finite dimensional distributions of the field Cy e + Z £ 
converge in law. Let k be some positive integer and xi,...,Xk points in M. d . For 
all £ = (£i, . . . ,£k) m (K d ) fc , we compute the characteristic function of (Cy e (xi) + 

z e (xi), ..., cy e (x k ) + z e (xk))- 

Q e (Q _ ^;/ e <Ei=i&-(oy«(«,0+««(*i))). 
By conditioning on the field generated by the white noise dWi, we get: 

Ce (£) = £;( e iCEj=i&-y«(*i) e -|;(E*=ife--F«(*i-i/)) 2 e* riJf!(v) - aCe *) 

Now, using proposition (13.11) . we have: 

K n rC it 

J2 / frFnixj - y)Q^(dy) ^ £ / ^(s,- - y)Q^(^), 

fdtti- F *( x i - ?/)) V- 271 W - fdtti- F *(*i - y)) 2 Q 271 (dy), 

J 3=1 J 3=1 

from where: 

QJg) ->. C(£) = E(e iCE ?= 1 ^^- Ffl( ^" y)Q7l( ^ ) e"^ (E ?= 1 ^- Ffl( ^" ?/))2Q271(c(?/) ). 

Thus, by applying Levy's theorem, we conclude that the finite dimensional dis- 
tributions of the field Cy e + Z e converge in law to those of a field X whose finite 
dimensional distributions are given by: 

#( e *l£=i &•*(**)) =C(0- 

Suppose that I is a positive integer that satisfies the condition of the proposition. 
For all £ in M. d , we have: 



We derive (13.61) similarly. 



E( e iCX{x) ) = E(e~ iC ^■ FR( - y)Q11{dv) e~^ {CFR{y))2Q2ll{dy) ). 

a 21 



We derive expression (13.51) by computing -J E(e i ^- X ^ x >)\^ =0 using proposition 13.11 



qAI 
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3.3. Scaling of X. The purpose of this subsection is to show that the field (X(x)) xeR c 
satisfies the multifractal scaling relation (this is what propositions 13.51 and 13.61 below 
assert). 

We first state two preliminary lemmas we will use in the rest of the paper. 

Lemma 3.3. Let S be some real number such that ^ 5 < a and 8 =£ a — 1. There 
exists C = C(5) such that we have the following inequality for \h\ ^ R: 

sup / \F R (y) - F R (y - h)\-±-dy ^ it^cA^ 1 . (3.10) 

Proof. By homogenity, we suppose that R = 1 and for simplicity, we suppose d ^ 2. 
Since j4j ^ t^ts + 1 and the right hand side of (I3.10P increases with 8, we have to 
show that for S G [0, a[ and \h\ ^ 1: 

sup f \ Fl (y)~F 1 (y-h)\ r ^-—dy^C\h^ M . 
x£R d JR d \ x y\ 

Indeed, this would imply that for \h\ ^ 1: 

sup / \F 1 (y)-F 1 (y-h)\ r -^—dy^C\h\^ M + C\h\ 
xem. d Jm. d \ x V\* 

<: 2C7|/i|( Q - 5 ) A1 . 

There exists C such that for all y and h, we have: 

\ip{y -h)- <p(y)\ ^ C\h\ and ip(y) ^ Cl\ y ^ 2 . (3.11) 

We set 

I(x)= [ \F 1 (y)-F 1 (y-h)\-^--dy. 
jR d \x-y\ d 

Therefore we get 

K x ) <C\h\ f . * 1 8 dy (3.12) 

J\ y \^3 \y\ F-r 
, n f i y y-h | i 

+ ^ / \-r~n — TT ~~ 1 771 — 77 1 H^V 

J\y\^i\y\ d - a+l \y-h\ d - a + ll \x-y\ s 

^C\h\+C [ I , ,/ 4 _ 1 - , y ~ Vj, * ,, dy, (3.13) 

where we denote by C different constants. 
First case: 5 < a — 1. 
Plugging inequality 

i 1/ V — h | (c? — a + l)|/i| 



1 lyld— a+l ly ^ld-a+1 I I'll /l| ci_a + :l - A |y| d— a+1 

in (I3.13p . we get 

I{x) < C|/i| / , TTT — ^ — , ,. — . — l -^- x dy. 

J\ y \ ^ 3 \y - h\ d - a+1 A |j/|*-«+i |x - # 
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We have: 

1 1 



dv 
V \^a\y- h\ d - a+1 A \y\ d - a+1 \x - y\ 5 

r l l r 11 

^ /, , Ui _ h\d—a+l \ry. _ q .\S y I \qi\d— a+1 I™ _ ,,|<5 " 

7|j/|^3ly Ai l l x y| J\y\^3\y\ \ x y\ 

f 11, 

<2SU P / | ? ,|d-a+l| r _ ?y |^ ^ 

which concludes the proof. 
Second case: S > a — 1. 

By the change of variable y = \h\u and setting h = \h\e with |e| = 1, we get: 

y-h y , 1 



i v i < 3 |y-/i| d - a+1 M d - Q+1 |x - y| 5 



iA< jl \u-e\ d ~ a+1 \u\ d ~ a+ll \x/\h\ -u\ 5 

U\ %; |, t | I I II I / I I I 

< |/i| a -« sup f i , M ~/ +1 -— ^_ i— L_^. 



a 

Lemma 3.4. Le£ <5 be some real number such that ^ 5 < 2a — d. There exists 
C = C(S) such that we have the following inequality for \h\ ^ R: 

sup / \F R (y) - F R (y - h)\ 2 -±-dy ^ C7|£|*~K (3.14) 

Proof. As in the proof above, we can replace |.|* by |.| and suppose that R = 1; thus 
we have to show inequality (13. 14j) with J(x) where we set: 

J(x)= f \F R (y)-F R (y-h)\ 2 -^—dy. 

Jut \x-y\ 

Using inequality (13.111) . we get 

J{X) < C ^ [ \vW-°)\x-v\' dy 

j\ y \ ^3 \y\ \ x y\ 
+ c f I y-h y , 2 i d 

J\v\<3 \y-h\ d - a+1 \y\ d ~ a+l1 \x-y\ s V 
i2 , „ f i y-h y l2 1 



^C\h\ 2 + C \, — *,. TT - T-rr -r h Tidy. (3.15) 

l\ y \^3 \y-h\ d - a+1 \y\ d ~ a+l1 \x-y\ s y v ; 
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Since 2 > 2a — d — 5, we only have to consider the second term in inequality (13.151) . 
By the change of variable y = \h\u and setting h = \h\e with |e| = 1, we get: 

]LzA y_A 1 dv 

\y\^\y-h\ d ' a+1 \y\ d - a+l1 \x-y\ s 

= \ h \^l-S f I U ~ e ?_* 1 du 

J\«\<a. \u-e\ d - a+1 \u\ d ~ a+l1 \x/\h\ -u\ 5 

< \h\ 2a - d - s snp [ | , U -\ -—^—f-^—du. 
l R d*\u-e\ d - a + 1 |w|^ a+l1 \a-u\ 5 



D 

Proposition 3.5. (Scaling along the even integers) 

Let I be an integer such that one of the following conditions hold: 

(1) I is even and l^\oJd < a — d/2 

(2) I is odd and (I + l)jlud < a — d/2. 

Let e be a unit vector (\e\ = 1). Then there exists Cf(e) > such that the following 
scaling relation holds: 

E{(X*{x + Ae) - W(x)f) ~ Cj{e){h««, (3.16) 

A— >0 it 

where we have 

( 2l = l(2a-d)-2j 2 1 u d l(l-l). (3.17) 

Proof. For simplicity, we will suppose that I is even and that IjfoUd < a — d/2. We 
introduce the following notation: 

f h (y)=Fj i (y)-Fi i (y-h). 

We shall see that the scaling at small scale of the sum (13. 6p is given by the term 
k = 0. Indeed for all k ^ 1 let us consider the integral 

fh(yi)---fh(y2k)(fh(y2k+i)) 2 ---(fh(yk+i)) 2 J J — : _, 2 „ ■ 



\ Vi-Vj \1l u d 
1 ^i<j ^2fc |— r— I* 



3 2 7 ?0(Q) P 4 7 ?0(^i) 



n 1'^,^ n l-m^> ^--- d ^- 

1 < i < 2fc I ^f - I * 2fc+l < i<j < fc+Z I ~r~ I * 

j>2fe 

< hJ k> h (3-18) 

where we set 



h= sup / |/fc(2/i)|...|A(3/2fc)| II 



j>2k 



2 X 

\ Vi-yj | 2 7i^d 
I u I* 



2ifVi-yj\ 

n e 7i<A(— FT 1 ) 
2— dyi . . . ci?/2fc- 
| 3/i-i/j |7i^d 
l^i<j <^2k \ R \* 
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and 
Jk,i= (A(j/2fc+i)) 2 • • • ifhiVk+i)) 2 - - .o dy 2 k+i ■ ■ ■ dy k+l . 

J{md\l-k , _?■■*■ ^, \ Vi—Vi pTi^d 

J \ K > 2fc+l s^i<j <fe+/ | R \* 

By using the estimates fl3.10n . fl3.14p and the inequalities ( 12.71) . (12. 8ft . one can show 
that for all fc ^ 1, we have 

I k J ktl < CR dk \^, 

with 

c M = (a-2(/-A;)7 1 2 ^)Al+((a-(2/-A;)7 1 2 ^)Al)(2A;-l) + (2a-d)(/-A;)-27 1 2 ^(/-A;)(/-A;-l). 

If a - 2(7 - k)>yjuj d < 1, then c k ,i = (2/ + fc(d - if^d); If a - 2(7 - k)^\uj d > 1 
and a — (2/ — k^\uJd < 1, then c^/ = C2Z + 1 — « + ^^ + (2/ — ?>k)^\ujd, otherwise 
Cfc^ = 2 A; + (2a — d){l — k) — 2^\u d {l — k)(l — k — 1). In all cases, it is easy to show 
that Cfc^ > (21 under the conditions of the proposition. 

Finally we study the term where k = 0. We get for h = \e with |e| = 1: 

4 7l ^(J^n) 



/ (A(yi)) 8 ...(A(w)) a n jxr^i---^- 



1 ^ i<n ^ / I jj 



w=Auj it J(Rd)' -^ «i — e| a a+1 R \ui\ a Q+i 



(<p(— (ui — e))i — j-3 — — Q?(— Mz)i — t-t — rr) 2 t? — du\...dui 

KYK R K "\Ui-e\ d - a+1 V R \uA d - a+lJ 1J - | A( Ui -u„) |4o>?a; d 

1 sj i<n ^ / I /j I* 



A-o K R J J {Rd)l K \ Ul -e\ d - a + 1 \ Ul \ d - a+l) '" K \u l -e\ d ' a+1 \ui\ d ~ a + lJ 



n- —5 — du\ . . . dur, 

1 < i<n < i I l n ' 

and inequality (12. 7p shows that this integral is finite when Ijfud < a — d/2. 

□ 
In the next proposition, we state the scaling relations of <Y along the odd integers. 
We define //(e) by: 

T j (p\= f ( u i~ eJ !^^_^2 / u t — e> u { 2 

* W " /, ,,/ li;, — P \d-a+l U.U-a+X' " ' M 7/ , _ p \d-a+l \„.\d-a+l> 

J (R d y \ a i e l l"i| |"i e i i u z| 



—5 — dui . . . dw/. 






Proposition 3.6. (Scaling along the odd integers) Let I be an integer satisfying the 
conditions in proposition ^. 5\ and 1 + 2l r yfuJd < a. 

Let e be a unit vector (\e\ = 1). Then we have the following scaling relation: 
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E{{X\x + Xe) - X 3 (x)) 2l+1 ) ~ Ef (e)(A)5H-i, (3.19) 

A— >0 it 

where we have 

( 2l+1 = l(2a - d) - 2 7 jW(Z - 1) + 1, (3.20) 

and 

S/(e)=7oC(/,7i)//(e)e J , C(Z, 7l )>0. 



Proof. As in proposition 13.21 setting C = 7o ° oe pd/2 -, it is possible to show that: 



^Coe-^iCi 



R d/2 
\2l+l\ 



E({X 3 {x + h)-X 3 (x)r +1 ) 

= V5 M C7 2fc+1 / f h ( Vl ) . . . f h (y 2 k+i)(fh(y2k + 2)) 2 ■ ■ ■ {f h {yk+i+i)) 2 x 

k=0 JiR^k+l+i 

11 ]yi -y 3 n^ d 11 | W - W |2 7l ^d 11 7^M^ dVl ■ ■ ■ dVk + l+1 > 

1 < i<j < 2fe+l | — ~ft~ \* 1 < i < 2fc+l | — ^f - | * 2k+2 < i<j < fc+Z+1 | — ^ — | * 

j>2fc+l 

where, as usual, we set: 

f h {y) = Fi{v)-F} l {y-h). 



Similarly to proposition l3.5l to get the main contribution as \h\ goes to 0, we examine 
the term k = 0. We introduce X: 



*rM^). 



x= / A(i/i)(A(i/ 2 )) 2 ...(A(i/m)) 2 n e "" ^ _ x 

y(R««)'+i /^ 2 

2^i<j </+l I R I* 

Putting h = Xe (|e| = 1), yi = i?«i, y, = Xu t (i ^ 2), we get: 

j „ i^/2 e aG-D^*(o)(A)S.H-i/i( e ) e . f V^(y)e 2l ^ p ^dy, 

where ^(l/) = l(l/) j i/ a+ i - Now a calculation gives: 

„2,.,_ /-oo 



2l-f(u; d f°° ! 2 . ? . w dp 



V^{y)e 2 ^^dy = -^ip( / r a "V(r)e 2 ^M ^dr)e,-, 

a Jq ar 

the last integral being negative since p(r) is a strictly decreasing function on the 
interval ]0, 2[: the result follows. Notice that the condition 1 + 2Z7 2 < a implies that 
this integral is finite while the other conditions imply that //(e) is finite. 

□ 
From proposition 13.61 we deduce readily that for A small the law of X{x + Xe) — 
X(x) is nonsymmetrical (for 7q 7^ 0). Indeed, by isotropy, we have : 

(X(x + Xe) - X(x)).e = X j (x + Ae,-) - X j (x) 

law 
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andSf(e i )=7 *Ca 7l )/f( e ,)>0. 

3.4. Tightness of X e and regularity of X. In this section, we prove that the 
convergence in law of X e towards X given by proposition 13.21 holds in a functional 
sense and that the field X is locally Holderian. The straight way to do so is to prove 
the tightness of the sequence X e by means of a Kolmogorov estimate (cf. chapter 13 
of [13]). 

Proposition 3.7. (Tightness) Let I be some positive integer that satisfies the con- 
dition of proposition \3.5\ and 7 a positive parameter such that jf < 7 2 . Then there 
exists e > and C independent of e such that for e < e and \h\ ^ R: 

Vx, E((X e (x + h)- X e (x)f) < ^pa-d^V^-i^ (3,21) 

and 

E((X e (0)f) <: C. (3.22) 

Proof. We only prove (I3.2ip (the proof of 13.221 is similar). We are going to compute 
the moment 

E{{X>{x + h)- Xl{x)f l ) = E(( [ U{y)^ {y) - C 'dW Q {y)f l ) 

where we set: 

f , s _ VR(y)y j vrjv - h )(y j - hi ) 

Je,h[y) \„.\d-a+l L. _ h\d~a+l 

\y\e \y n \e 

We get: 

E(( [ Uk(y)e xt(y) - c <dW (y)) 21 ) = e" 2 ^ / f e , h ( yi ) . . . f e , h {y 2l )E{e* € dW { yi ) ■ ■ ■ dW (y 2l )), 

(3.23) 
where 

x* = x e { yi ) + ... + x*{ y2l ). 

The rest of the computation can be performed rigorously by regularizing the white 
noise dWo, using lemma I2T21 and going to the limit. It is easy to see that we obtain 
the same result by introducing the following formal rules: 

E(dW {y)dW (y')) = 6 y _ y ,dy (3.24) 

and 

E(dW (y)X e (y')) = l0 (e)k?(y' - y)dy (3.25) 

As a consequence of lemma 12. 2} E(e Xi dW (yi) . . . dW (y 2 [)) is the sum of terms of 
the form 

E(dW ( yi )X e ) . . . E(dW (y k )X*)E(dW (y k+1 )dW (y k+2 )) . . . E(dW (y q . 1 )dW (y 2l ))e^ ( ^ 2 \ 

(3.26) 
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We will compute the limit of each one of these terms. By using (13.251) . we get 

2/ 



E(dW (y k )X*) = 7 o(e)(]T fcf ( Vi - y k ))dy k 
= lo (e)k?(0)(l + Qi)dy k 

where 

« = ^o)(E*f(w-y*))- 

We also have from the definition of X e : 

e §£((* e ) 2 ) = e ap £/fl (0)+E l < J P e /i i (^ 1 ))(70(e) 2 +7 1 2 ) > 

By using lemma |2~72| expression (13.231) and the rules above, we get: 
£((/ H ,/^(y)e^^<W (y)) B ) = ^ 



Imd) k +i fe,h{yi) ■ ■ ■ Uh(y2k)(fe,h(y2k+i)) 2 ■ ■ ■ (/ e ,h(i/ft+j)) 2 n^il 1 + QU,i) e k '' d vi ■■■ d v, 



k+l 



where 



< 2A; j>2fc 



and 



^ = ((7o( e )) 2 +7 1 2 )( E p^(^)+2 E E^(^: 



1 < i<3 <, 2A; 1 < t < 2fc j>2k 

+4 £ ft/M^))- 

2fc+l < i<j ^ fc+/ 

We first take care of the normalizing constant outside each integral: 

(7 (e)A;f(0)e" 1/2((7o(e))2+7l2)ft / fl(0) ) 2fc e 2Z((7o(e))2+7l2)Pe / fl(0) " 2 ' c,e . 

By the choice of C e , we have e^iii^f+iDp./R^^ic, = 1 Using expansions ([23]) and 
(12. 5p . we derive the following limit: 

Cne- 1 / 27 ^ 1 



7o(e)^(0)e- 1/2((7o(e))2+72) ^/« (0) -> 7o 



In conclusion, the constant outside the integral of term k in the above sum converges 

to«, ,r 7 " Coe " 1/27?G M 2fc 

ro a £:A r372 J • 

Let 7 be such that jf < 7 2 . One can choose eo > such that 7o(eo)) 2 + 7i < 7 2 - 
Using the fact that, for all y, p e /n(y/R) ^ uJd^i + -^ + C with C independent of e, 
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we get: 



e su <C TT L_ TT 1 , 



3/i-Vj 17 ^d 11 I j/j-j/j |27 2 o; d 
1 < i<j < 2fc I i?, I* 1 < i ^ 2fc I R I* 

j>2fc 

is— ■ (3-27) 

i y»-yj |47 2 ^d v y 

2/c+l ^ i<i <fc+z I i? I* 

Finally, we conclude by using the fact that |Qffcj| is bounded by a constant inde- 
pendent of e, inequality ( 13.101) and ( 13.141) similarly as in the proof of proposition 

D 

Corollary 3.8. One can easily deduce from this that for jf sufficiently small, by 
Kolmogorov's compacity theorem (\13\), X t tends to X in the functional sense and 
that X is locally Holderian. 

Comment 3.9. Starting with a two parameter (R, a) monofractal Gaussian field, 
we constructed a four parameter (R, a, 71, Jq) multifractal field with nonsymmetrical 
increments. This family has it's own interest. As we shall see in the next section, 
this family is too restricted to take into account all the constraints needed for a 
satisfactory model of turbulent flows. 

In the case where Jq = 0, we obtain symmetrical random fields which extend to 
higher dimensions the model introduced in [2]. 

In the next section, we will study a multifractal field which is not in this family 
but that can be seen as a limit case where 71 = and 70 is constant (independent 
of e). As we will see, this family will be compatible with the 4/5-law. 

4. A STEP TOWARDS A MODEL OF THE VELOCITY FIELD OF TURBULENT FLOWS 

An acceptable solution to the problem of hydrodynamical turbulence in dimension 
3 would be to construct a random velocity field U solution to the dynamics (Euler 
or Navier Stokes typically) that is stationnary, incompressible, space-homogeneous, 
isotropic and that satisfies the main statistical properties of the velocity field of 
turbulent flows. Two main properties are: 

(1) The 4/5-law of Kolmogorov that links the energy dissipation of the turbulent 
flow to the statistics of the increments of the velocity. This law is widely 
accepted since it is the only one that can be proven with the dynamics ([6], 
[7], [II]). More precisely, this law states: 

E((U(x + 0-U(x).-^) 3 ) = ~D\Z\. (4.1) 

In the above formula, D denotes the average dissipation of the kinetic energy 
per unit mass in the fluid. 

Remark 4.1. To obtain this law, it is sufficient to suppose that the field U 
is space homogeneous and isotropic. 
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(2) The intermittency of the field U: 

E((U(x + - U(x).^r) ^ C q \^, (4.2) 

where ( q is a well known concave structure function (cf. [7]). 

It is a very challenging task to construct a field with all the aforementioned prop- 
erties, especially because this field must be invariant by the Euler or Navier-Stokes 
equation. 

Nevertheless, one can in the first place forget the invariance by the dynamics and 
simply try to construct a field that satisfies the other properties. The 4/5-law shows 
that the nonsymmetry of the increments is an essential feature. Let us consider the 
family X constructed in the previous section (d = 3). By proposition 13.61 we have: 

E{{{X{x + Ae) - X{x)).ef) ~ tf s (A)& C 3 > 0, 

A— >0 It 

with £ 3 = 2a — 2. To satisfy the 4/5 law one should have £3 = 1, which gives 
a = 3/2. This is incompatible with the constraint 3/2 < a < 5/2. Thus we have 
now to modify the family X to reach the limit case £3 = 1. In this aim, we will 
construct a new (three parameter) family X corresponding to the limit case 71 = 0, 
70 constant > 0. 

4.1. Construction of the field Xq. In this section, we only outline the main steps 
of the construction of Xq. The field Xq j£ is given by formula (I2.3P where X e is now 
defined by: 

X € (y) = lo [ kf(y - a)dW (a). 



We suppose that a is in the interval ]0, 1[. We choose the normalizing constant C e 
such that: 

7ofcf (0)e- C£+ ^oft/fl(o) = 1. 

We start by stating a lemma we will use in the proof of the proposition below: 

Lemma 4.2. let 5 be some real number different from d. Then there exists C = 
C{5)>0 with: 

d± ^ Ce (d-6)A0_ (43) 

'\u\ ^R Pie 

Proof. We suppose 5 > d, the other case being obvious. We have: 

du d _ s f du 



»KR \u\e «=™ J\u\^R/e {f\ v] ^i6(v)\v + u\dvY 



\s 

"- Cli J\u\ ^R/e VJ 

du 



<e d - s 



K d (f\ v \^i0( v )\ v + u\dv) 



D 



We can now state the following proposition: 
Proposition 4.3. Let I be an integer ^ 1 and 70 such that: 
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(1) J$U d < a if 1 = 1. 

(2) (21 - 3/2)7^ <aAd/2ifl>l. 

Then for all x, X 0;6 (x) converges in L 21 to a random vector X (x) (i.e. E((X 0j6 (x) — 
Xo(x)) 21 ) — > 0). The random vector field Xq(x) satisfies the following scaling: For e 
(\e\ = 1), and for all q ^ 21: 

E{{XJ{x + Ae) - X>(x)Y) xZq Cj(e)(A)^, (4.4) 

where ( g = qa — \q(q — l)7o^d an d 

C j( e) = R qd/2 Skr21 7 2 m J" -pr 1 -pr u\ u\ - e> 

(4.5) 



E{{xi{x)T) - / n ^? TT r#¥ ^i---^- (4.6) 



Proof. We will first prove that: 

eggg -pr MVi)yl 

\ Vi-Vi |To^d J-J- 
1 < *<i < 2/ I ji I* l^i<2Z 

We remind that the right hand side of the above limit exists by lemma 12.41 In order 
to prove the above relation, we develop E((Xq e (x)) 21 ) in I + 1 terms similarly as 
in the proof of proposition 13.71 then, using formula (j2.4p and the fact that, for all 
y, p e /n(y) ^ cu^ln — + C, we are led to show that, for all k ^ I — 1, we have the 
following convergence: 

e (/-fc)(d- 7( > d ) e -2(«-fc)(i-fe-l) 7o 2 a; de -4fc(«-fe)7 2 ^ f 'PrJvM fRJ^Mk „ 

•W+. |yi|f- a+1 '" li/2fclf- Q+1 

(^(2/2fc+i)^ fe+1 ) 2 {(pR{yk+i)yi+i) 2 i-r e^^V 1 ) 

— ,2(^+1) ••• , ,2(^+1) 11 — ^%---^^0' 

|y2fc+l|e |i/fc+«|e 1 <«<j ^2fc I ji I* 

We apply inequality (14.31) and obtain (if a = d/2, one can work with a — r\ for rj > 
sufficiently small) : 



JR d \v\ e 



iR d \y 

Therefore the above convergence to amounts to showing that, for all k ^ I — 1, we 
have the following inequality: 

d + (2a - d) A - -iluo d > 2(1 - k - l)7 2 cj d + Ak^ooa. 

This is equivalent to (21 — ^) r y 2 uj d < a A |. One can show, for all x, that (A' 0e (a;)) e> o 
is a Cauchy sequence in L 21 by computing E((X^ e (x) — Xq € ,(x)) 21 ) and letting e,e' 
go to 0. Thus, E((Xq(x)) 21 ) is given by the left hand side of (14. 6ft . 
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To show the scaling (14.41) . observe that we can prove the following analogue to 
flUD for any q ^ 21: 

r _ eiM^) r- 

E{{X>(x + Xe)-X>(x))*)= II — - — a- J] he(yi)dyi...dy, 

• / CR d )« 1<i<J - <g |^|Z d i<^ 9 



(4.7) 



where 

^fl(y)y J _ <?fl(y- Ae)(y J - Ae J ) 

|z/| d_a+1 |l/ - Ae| c 

By setting y; t = Xui in the integral of (14.71) . we deduce easily (14.41) . D 



f (,.\ _ ^R\y)y __ ^R\y ~ ^ji;/ - ^ c ; / , < 

Ae^J — I..M_a+l |„ _ \ P \d-a+l ' ^ 



Remark 4.4. It is not obvious why in the above proposition the coefficients CUe) 
are different from (cf. appendix). 

Remark 4.5. Similarly as in the previous section, for 70 sufficiently small, Xq^ 
converges in law to X in the space of continuous fields. 

4.2. Nonsymmetry of the increments of Xq. By isotropy, we have: 

(X (x + Ae) -X (x)).e = X$(x + Xej) - X$(x) 

law 

One can show that for A small the law is nonsymmetrical by showing that the 
third moment is 7^ 0, that is C 3 z {ef) 7^ (see appendix). 

4.3. Towards a model of the turbulent velocity field. In dimension 3, for the 
field X , we have: 

( q = qa- 2-nq(q - l)^. 

Thus, for a = 1/3 + 47T7q, we have £3 = 1, which means that for this choice the asso- 
ciated fields X satisfy at small scale the 4/5 law with a non zero finite dissipation. 
Unfortunately, the fields in this family are not incompressible. The incompressible 
case (at small scale) would correspond to the choice a = 1, a limit case which is 
excluded by the constraint < a < 1 needed for the validity of the scaling of propo- 
sition H]3j There is another severe obstacle for the choice a — 1. Indeed, for the field 
Xq above, we have: 

C g = (l/3 + 67T7 2 )g-27r 7o V. 
One can easily identify the intermittency parameter 47T7q using the experimental 
curve given in [I] (cf. fig 8.8 p.132 in [7]). With this data, we find iirj^ = 0.023. 
With this small intermittency parameter, we would get a ~ 0.35 which is not close 
to the incompressible value a — 1. So, in spite of its qualitative interest, this model 
cannot reach quantitative adequacy. 

Another natural way to get incompressible fields is to use a Biot-Savart like for- 
mula and take the limit as e goes to of fields of the form: 

W{x) = [ ^ R {x-y) X ~ y _ a+1 AdQ% 
Jw.3 \x-y\i a+1 

where dQ t is an isotropic random field. For example, we can take: 

dQ e = e xciy) - a dW{y), 
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where dW(y) = (dWi(y) , dW 2 (y) , dW 3 (y)) denotes a three dimensional white noise 
and X e is defined by the following formula: 



X%)=7 / Kf{y-a).dW{a) 

JR3 



with K R (x) = j^jj2l\ x \^R. 

As for the case of X , we choose the constant C e such that U e converges to a 
non trivial field U as e goes to zero. The vector field U we obtain is incompressible, 
homogeneous, isotropic and intermittent with structural exponents ( q given by: 

( q = qa - 2vT7 2 g(g - 1). 

Unfortunately, since the field dfl e (y) is isotropic with respect to all unitary trans- 
formations (and not just the rotations) we get for U the symmetry: 

U(-x) - 17(0) = U(x) - £7(0) 

law 

so that the dissipation is equal to 0. Thus the construction of an homogeneous, 
isotropic, intermittent and incompressible vector field with positive finite dissipation 
remains an open question. 

Comment 4.6. In our approach, we perturb a Gaussian field to get multifractality 
and we further introduce some dependency to obtain also dissymmetry. This can 
make one think that in turbulence dissipation is linked to intermittency. This is a 
rather intricate issue. On one hand, only dissymmetry seems to be needed to get 
energy dissipation (see the 4/5 law). On the other hand, it is well known experimen- 
tally that dissipation is not homogeneously distributed in the fluid but rather follows 
the lognormal distribution described by Kolmogorov and Obukhov by which it appears 
linked to intermittency. 

In dimension d = 1, our model displays some kind of (non causal) leverage effect. 
To get a realistic model for finance, with causal leverage effect, we have to make some 
specific change in the construction. This issue will be addressed in a forthcoming 
paper. 

5. APPENDIX 

In this appendix, we prove that for q even C*(e), given by equation (14.51) . is 
different from outside a countable set and in the neighbourhood of 0. We also 
show the same result for C^(ej). 

Consider first the case q even; we set q = 21 with I greater or equal to 1 and we 
introduce the following function F: 

f(i)= f n i — 1 —^ n f(ui)d U1 ...du 2h 

where / is the real function defined by: 

u j + e j /2 u j - e j /2 

* ^ ~ \u + e/2\ d ~ a + 1 ~ \u-e/2\ d ~ a + 1 ' 
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The function F is analytical in a neighbourhood of 0; therefore, in order to obtain 
the desired result, we have to prove that F is not identically equal to 0. One can 
show that, for all i < I, F®(0) = and that: 

o/| r i 

F(0 (°) = ^( / ln (i ,)/M/M^icM'- 



The Fourier transform of Mm) is adPf(jFTd) + bdS where a^ > and bj are two 
constants that depend only on the dimension and Pf is Hadamard's finite part (see 
p. 258 in [15j). Since J Rd f(u)du = 0, we get: 

1 f /(f) 2 

Ml r)f{ui)f{U2)duidU2 = a d -J7UT d ^ 

thus F«(0) > 0. 

Let us now consider Cg(e_j) and the corresponding function: 

^(7) = / I m iri -f(u 1 )f(u 2 )f(u 3 )du 1 du 2 du 3 . 

\3 \u\ - u 2 V 1 \u\ - u 3 P \u 2 - W3 7 



We obviously have F(0) = 0, F'(0) = and: 

-F"(0) = 1=1 ln(|«i - u 2 |) ln(|wi - u 3 \)f{ui)f{u 2 )f{u 3 )du 1 du 2 du 3 

so that X = J" Rd f(x)Q(x) 2 dx where: 

e(x) = / Mk-i/D/d/)^- 

Now we prove that there exists some real constant c different from such that: 

8 (a; = c (1 1 —7- — . 5.1 

v ; y \x + e j /2\ l ~ a \x-e j /2\ 1 ~ aJ v ; 

Indeed, we have (in what follows, c denotes different real constants that are not 
equal to 0): 

M) 



9(0 = c- 

"I" nn c 



and /(£) = csin(7r^) F] |+ T thus 



e(0 = csin(7rf 



The above expression (15 .ip now follows from: 

xi £i 



and 2isin(7r^) = <L e ./ 2 (f) - &<% /*(£)• 
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Now let us denote x = x^ej + x and 

., , S/ + 1/2 2/~ 1/2 



((y + 1/2)2 + a)P ((y _ 1/2)2 + o)? 

If we set: 

/ i\ . ,d — a + I ■ ._ l2 . 
¥W = 0( ^ ,a^, Fl ) 

and 

^<y) = 0( — ^ — ,^ j , is"i 2 ), 

we get: 

I = c 2 dx / (p(x^)i(^(x^)' 2 dx-' 



Since < a < 1, it is easy to check that V^- 2 ) is a positive function of z, decreasing 
on [0, oof. One can also check that there exists some z* > 1/2 such that tp(z) is 
positive on [0, z*[ and negative on ]z*, oof. Since <p and ip are even and f Q <p(z)dz = 0, 
one can derive the following: 

poo 

(p(z)^(z) 2 dz = 2 ip(z)^(z) 2 dz 
Jo 

pz* />oo 

= 2/ ^(z)^(z) 2 dz + 2 <f(z)ip(z) 2 dz 

J0 Jz* 

fZ* /"OO 

^2 ip(z)ip(z) 2 dz + 2^(z*) 2 <p(z)dz 

Jo Jz* 

= 2 [ Z ip(zmz) 2 - 4>(z*) 2 )dz 

Jo 
>0. 



It follows that F"(0) > 0. 



□ 
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